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The unequivocal detection of entanglement between two distinct matter-wave pulses is a significant challenge 
that has yet to be experimentally demonstrated. We describe a realistic scheme to generate and detect continu- 
ous variable entanglement between two atomic matter-wave pulses produced via degenerate four-wave-mixing 
from an initially trapped Bose-Einstein condensate loaded into a one-dimensional optical lattice. We perform 
a comprehensive numerical investigation for fixed condensate parameters to determine the maximum violation 
of separability and Einstein-Podolsky-Rosen inequalities for field quadrature entanglement, and describe and 
simulate an experimental scheme for measuring the necessary quadratures. 
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I. INTRODUCTION 

Recently there has been much interest in the quantum prop- 
erties of matter waves and the study of quantum atom op- 
tics illl^l- The field has grown out of atom optics and encom- 
passes concepts and ideas from quantum optics, condensed 
matter theory, atomic and molecular physics and quantum in- 
formation theory. Experimentally, ultra-cold atoms provide a 
clean and controllable environment to investigate a wide range 
of new and existing models. Notable examples include the ob- 
servation of the BCS-BEC crossover regime jsl and the Mott 
insulator-superfluid quantum phase transition |4|] . 

In addition to answering fundamental questions of science, 
this developing area of physics has the potential for new and 
exciting applications. It is predicted that quantum entangle- 
ment will enable a novel set of technologies based on fun- 
damental quantum principles, such as precision measurement 
devices [5;] and quantum computers la]. The precise coher- 
ent quantum manipulation of ultra-cold atomic systems has 
been demonstrated in many experiments and holds promise 
for these systems to be candidates for quantum information 
applications in the future. 

There have been a number of demonstrations of quantum 
atom optical phenomena in recent years. These include the 
observation of non-classical effects in atomic fields such as the 
Hanbury Brown-Twiss effect for bosons anti-bunching 
for fermions i^, reduced local pair correlations ifiol [Till . 
sub-Poissonian number fluctuations lfT2l [TbIi. a nd density cor- 
relations from molecular disassociation lll4ll . atomic colli- 
sions [[ist], and in the Mott-insulator regime in an optical lat- 
tice 11 1611 . The majority of these investigations have concen- 
trated on coiTelations in the atom number or density. Future 
applications of quantum atom optics utilizing entanglement 
and squeezing will require manipulation and detection of the 
phase of the quantum state of matter-waves. This presents a 
challenge as creating stable phase references and performing 
mode-matched intereference is likely to be difficult for ultra- 
cold atoms. 

Entanglement can be generated and utilized in a variety of 
forms. In this paper we address the generation and detection 



of continuous-variable entanglement between orthogonal spa- 
tial modes of an atomic field in a second-quantized formalism. 
This is distinct from entanglement between single particles, as 
seen in a first-quantized picture. Such continuous-variable en- 
tanglement has been extensively used in the field of quantum 
optics IitIi . 

Entangled and squeezed states of the electromagnetic field 
have been generated in quantum optics experiments, and 
proof-of-principle demonstrations for potential applications 
such as quantum cryptography have been successful IitIi . 
Photons may not be ideal for all quantum operations; light 
travels fast but is difficult to contain, and the lack of mass 
makes it relatively insensitive to rotation and acceleration in 
interferometric sensors |18-]. In principle, the achievements of 
quantum optics can be replicated with bosonic atoms. Indeed, 
there have been impressive achievements in entangling the 
spins of two distinct atomic clouds in the continuous-variable 
limit ifTgll . Recently, entanglement in the reduced two-body 
density matrix was demonstrated by number and phase cor- 
relations in double- and few-well systems ifisll . However, no 
experiment has yet demonstrated entanglement between the 
spatial modes of an atomic field. 

Several suggestions have emerged for creating entangle- 
ment in a Bose-Einstein condensate (BEC). One possibility 
is to transfer entanglement from an pre-existing entangled 
source to the atoms. For instance, the state of entangled 
light can be mapped onto atomic modes by using lasers to 
induce a Raman transition i20ll2lll . A related suggestion for 
entanglement detection involves mapping the atom statistics 
onto photons in the reverse process, in order to access well- 
developed sing le photon and quadrature measurement tech- 
niques ll2lll22ll . 

A second approach is to generate the entanglement directly 
in the atomic system. Atomic systems exhibit three- and 
four-wave mixing processes that can generate entanglement in 
light. The dissociation of a diatomic molecular Bose-Einstein 
condensates into entangled pairs is analogous to the process of 
three-wave mixing in an optical parametric amplifier ll23il24ll . 
In the second-quantized picture the spatial modes of the prod- 
ucts become entangled. Unfortunately, the entanglement that 
results is difficult to use or detect because of the necessity to 
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use a phase reference that is correlated to the phase of both 
the atomic and molecular modes. Atom-light scattering is an 
example of a four-wave mi xing process which can entangle 
atomic and photonic modes llzsll . Superradiant Rayleigh scat- 
tering occurs in an elongated condensate, where the scatter- 
ing process is self-stimulating, generating only one or a few 
macroscopically occupied, entangled modes ll26l l27ll28ll . 

The entanglement we consider in this paper occurs between 
modes of the atomic field. The inherent s-wave scattering be- 
tween ultra-cold atoms can result in the four-wave mixing of 
matter waves. Stimulated four-wave mixing between colliding 
condensates was first reported in 1999 12911 , followed by ob- 
servations of spontaneous four-wave mixing and the demon- 
stration of the potential reversibility of the process ll30ll . In 
principle, this four-wave mixing scheme produces many pairs 
of spatially entangled modes Oil l32ll . in analogy to a recent 
optical experiment [33t|. However, the sheer number of res- 
onant modes limits the gain and entanglement generated be- 
tween any given pair of modes, as the condensate is rapidly 
depleted. In this work we extend a suggestion to use degen- 
erate four-wave mixing (where the two input waves are the 
same mode) of a quasi-lD BEC in an optical lattice ll34i Issll 
to create and detect continuous variable entanglement. The 
quasi- ID geometry limits the number of modes for resonant 
collisions, and leads to the generation two separate yet entan- 
gled matter waves. 

Hilligs0e and M0lmer Js^l first suggested loading a station- 
ary Bose-Einstein condensate into a moving optical lattice in 
order to make a degenerate four-wave-mixing process reso- 
nant. There is a large body of work on the inherent dynamical 
instabilities in moving optical lattices which cause heating of 
BECs (see ll36ll37i ISSIl and references therein). Interestingly, 
the same dynamical instability can be interpreted as a coUi- 
sional process that generates entangled modes of specific mo- 
menta. An experiment was subsequently performed by Camp- 
bell et al. [35,] that shows pairs of modes are populated with 
the momenta predicted by mean-field theory; however there 
has been no experimental proof of entanglement in this sys- 
tem. The goal of this work is to propose a method to demon- 
strate entanglement between modes populated by degenerate 
four- wave-mixing . 

To utilize or detect the presence of entanglement requires a 
suitable phase reference lll7l[39ll40ll which is provided by the 
near classical coherent output of a laser in optics. In atom op- 
tics the equivalent of laser light is a Bose-Einstein condensate. 
However, in experiments the size of BECs is typically limited 
to be between 10"^ and 10^ atoms, and interactions between 
atoms result in atom losses, number dependent phase evolu- 
tion, and phase diffusion meaning that they are less than ideal 
as a phase reference. Also, setting aside a separate conden- 
sate of atoms to use as a phase reference may not be practical 
experimentally. 

In Ref. ||4lll we extended the separability criterion of Duan 
et al. ll42ll and Simon ||43}l , and the criterion for demonstration 
of the Einstein-Podolsky-Rosen (EPR) paradox of Reid lliill 
to situations where no assumptions could be made about the 
quantum state of the available phase reference. This allows 
the possibility of using a non-classical local oscillator in bal- 



anced homodyne measurement schemes for quadrature mea- 
surements that are common in quantum optics. We then ap- 
plied these criteria to a toy four mode-model of degenerate 
four-wave mixing in a Bose-Einstein condensate with up to 
1000 atoms, and showed that appropriate beam-splitting oper- 
ations combined with atom counting could be used to demon- 
strate both inseparability and the EPR paradox. 

In this paper we extend this approach to perform one- 
dimensional simulations of a more realistic experimental sys- 
tem. We simulate the adiabatic loading of a trapped Bose- 
Einstein condensate into a moving optical lattice to initiate 
degenerate four-wave mixing, before turning off the lattice po- 
tential and performing beam-splitting operations with appro- 
priate Bragg pulses to generate four distinct atomic clouds. 
Performing number difference measurements between these 
clouds provides access to the atomic quadratures and the sub- 
sequent violation of the entanglement inequalities. We de- 
termine how the measure of entanglement depends on the 
time held in the optical lattice, the relative phase of the 
Bragg pulses, and the number of seed atoms in the outgoing 
modes. We conclude that these experimental complications 
will not necessarily prevent the demonstration of entangle- 
ment of matter-waves in such a system. 

The paper is organized as follows. In Sec.|II]we review the 
degenerate four-wave mixing of a BEC in an optical lattice 
that was the subject of previous theoretical and experimental 
investigations ||34[ [35l l45ll . Section |III] details the proposed 
measurement scheme, including the basics of Bragg pulses, 
homodyne measurements and the entanglement criteria we 
employ. The results of numerically implementing our scheme 
are presented in Sec. llVlbefore we conclude in SecIV] 



II. DEGENERATE FOUR- WAVE MIXING OF 
MATTER-WAVES 

The dominant interaction between ultra-cold atoms in a 
Bose-Einstein condensate is s-wave scattering. In free space 
the conservation of energy and momentum in the collision 
of two particles requires that the outgoing particles have mo- 
menta falling on two opposing points of a sphere in fc-space in 
the centre of mass frame. Such s-wave scattering spheres have 
been observed by colliding two BECs ifisilsoll . Alternatively, 
this can be viewed as a four-wave mixing process which is 
resonant and phase-matched on the s-wave scattering sphere. 

The presence of a periodic potential alters the dispersion re- 
lation for the particles and thus the resonance conditions for 
collisions. Hilligs0e and M0lmer ll34ll showed that adiabat- 
ically loading a BEC into a moving optical lattice can allow 
degenerate four-wave mixing, where the two incident particles 
are in the same initial state. The optical lattice is created by 
shining two lasers with wave-vectors ki and k2 of similar fre- 
quencies LOi and u}2 onto the atomic cloud. The light creates 
an effective potential for the atoms 

y(r) = -ysin(2ki-r-,5t + 0), (1) 
where 2ki = ki — k2, (5 = — L02, is the relative phase 
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FIG. 1: (Color online) Band structure of an optical lattice with s — 
1. Two atoms in the (moving) condensate mode can collide into 
modes 1 and 2, conserving energy and quasi-momentum. 



between the beams, and Vl is proportional to the intensity of 
the lasers and depends on the details of the atom-light interac- 
tion. It is convenient to write Vl = sE^^ where s characterizes 
the strength of the optical lattice, and Ejj = /i^fc|/2m is the 
recoil energy. 

The degenerate four-wave-mixing process is illustrated in 
Fig.[T]where two particles in mode with quasimomentum qq 
in the lowest Bloch band can colUde, transferring an atom to 
each of modes 1 and 2 with quasimomenta qi and (72 ■ During 
this process they conserve their total energy 2eo = £1 + £2 
and quasimomentum 250 = qi + q2 modulo 2fci. By subse- 
quently adiabatically turning off the optical lattice the gen- 
erated quasi-momentum states are converted to momentum 
states in free space between ifci. If we write resulting mo- 
mentum of mode I as fci = fco — Afc, then the momentum of 
mode 2 is ^2 = fco + — 2kL. 

Defining the annihilation operator of mode j to be the 
model Hamiltonian describing this process can be written as 

H = ihx (jiQa\a2 — 0,^^0,-^^0,2^ ■ (2) 

This pairwise scattering process is predicted to produce sub- 
Poissonian number correlations and quadrature entaiiglement 
between the atoms with quasimomenta qi and q2 ifsil i45ll . 
However, the measurement of number correlations or demon- 
stration of entanglement was not accomplished in the sole ex- 
perimental realisation of degenerate four-wave mixing in a 
EEC to date Issll . 

We can perform a simplified analysis of the Hamiltonian 
in Eq. (|2]i using the undepleted pump approximation where 
we assume that mode begins in a large coherent state 
where a is real and the number of particles A^o = ct^- At 
small enough times one can assume the population of mode 
is undepleted and therefore the Hamiltonian can be approxi- 
mated as 

H w ihNox (^o,\a\ — 0^0,2^ ■ (3) 
In the Heisenberg picture the solution for this system is IitI 



a-^{t) ^ cosh{Noxt)a-^iO) + sinh(iVoxO«2(0), 
02(0 = cosh(7Voxi)a2(0) + smh{NQxt)a\{0). (4) 

If modes 1 and 2 are initially vacuum then they make up a 
two-mode squeezed state that is well-known from quantum 
optics ll46ll . The two modes are exactly correlated in number, 
and anti-correlated in phase. Measuring the correlation be- 
tween the number of atoms in modes 1 and 2 is, however, in- 
sufficient to prove entanglement, as it is possible to construct 
a separable density matrix that is consistent with any set of 
local number measurement outcomes ll40ll . 

To unequivocally demonstrate that entanglement exists be- 
tween the two modes it is necessary to perform phase sensi- 
tive measurements. The phase difference between the modes 
is insensitive to anti-correlations, and thus measurement of 
the relative phase will not help demonstrate entanglement be- 
tween these states. We therefore need to measure the phase 
quadratures of both modes. 

Quadrature measurements have been used to prove entan- 
glement between photonic modes in quantum optics experi- 
ments 03] • A third phase reference, or local oscillator, is in- 
terfered with the entangled modes on a beam splitter to reveal 
the quantum coiTelations. A large and coherent local oscillator 
(produced by a laser) allows one to make accurate measure- 
ments of the field quadratures, defined as 

X^=o^+a], Y^=i{a^-a]). (5) 

Duan et al. ||42ll and Simon [43] have derived a simple criteria 
for the separability of the two modes. All separable states 
obey 

Var [X^ - + Var [Yi + ^2] > 4, (6) 

where we write the variance Var [A] ~ {A^) — {A)"^. Vio- 
lation of the above inequality indicates that the system must 
be entangled. For all times i > the solution Eq. (|4|i of the 
Hamiltonian (|3]l in the undepleted pump approximation vio- 
lates this bound 

Var[li - X2] + Var[fi + Fj] = 4e-2^«*. (7) 

This inequality is useful for detecting the entanglement gen- 
erated by the four-wave mixing process, and has been suc- 
cessfully employed in optical experiments UtIi . The Einstein- 
Podolsky-Rosen paradox has been demonstrameted in optical 
fields using similar inequalities iflTi l44ll . 

III. ENTANGLEMENT DETECTION 

In optics the quadrature measurements required to demon- 
strate entanglement can be made by balanced homodyning 
that mixes the signal beams with a phase reference beam on 
a 50-50 beam splitter We employ balanced homodyning for 
atoms as this method typically has a superior signal-to-noise 
ratio compared to unbalanced homodyne schemes. In the 
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Heisenberg picture we define a = a^^ and b = b-^^^ as the 
annihilation operators for the signal and local oscillator be- 
fore the beam splitting. The output modes a^^t and b^^^^ are 
given by 

Oout = « + U/^, and foout = (flin - ^i„)/^- (8) 



The final step is to measure the difference in the number of 
particles exiting the beam splitter ports a^^j and b^^^^.. The 
measured quadrature is rescaled by the size of the local oscil- 
lator, according to 



X = 



b 

"out "out 



(6^ I )V2 
\^in in/ 



b)b + a6t 
(&t6)i/2 



+ a, 



(9) 



where the final approximation is good for phase references 
that are large coherent states. 

Thus three ingredients are required to perform homodyne 
measurements with Bose-Einstein condensates: a suitable 
phase reference, the atomic equivalent of a beam splitter, and 
detectors able to accurately count the number or density of 
atoms in the various modes. Bragg scattering can be used to 
interfere (beam split) atoms of different momenta, and rea- 
sonably efficient atom detection has been demonstrated with 
multi-channel plate detectors for metastable Helium ifTilgl fTsll . 
ionized ^^Rb detection ll47ll . high numerical aperture optical 
techniques 1480 . and utilizing optical cavities |8i]. Combined 
with time-of-flight expansion, the population distribution in 
momentum space can be directly measured. 



A. Bragg pulses 

Bragg scattering is the process of applying a moving opti- 
cal lattice to coherently transfer atomic populations from one 
momentum state to another The potential in Eq. ([T]) will reso- 
nantly transfer an atomic population with momentum ki to 
k2 (and vice-versa) if and only if ki — k2 = ±2ki and 
5 = fi-k^ • (ki + k2)/m; this is equivalent to assuming the 
atom scatters a photon from one beam to the other, conserv- 
ing both energy and momentum. Higher order scattering is 
possible, but will not occur if the lattice is relatively weak 
is < 1). 

A weak lattice held for a duration r = Trh/sEn provides 
a 7r/2 Bragg pulse, and is equivalent to a 50-50 beam splitter 
operation. This can be most easily seen in the Heisenberg 
picture 



exp(— ZTr/c^/fc^s) 



a2{t + t) — cxp(— jTr/cj/fc^s) 



^2 ft) - ai{t)e~ 
V2 



(10) 



.(11) 



where the pulse begins at time t. Note that the relative phase 
of the two laser beams 9 contributes to the evolution, and this 
is important when considering quadrature measurements. 

In this paper we will make use of mode as a local oscilla- 
tor in order to detect entanglement between modes 1 and 2. As 



simultaneous quadrature measurements of modes 1 and 2 are 
necessary, we require two local oscillators for balanced ho- 
modyning. This is achieved by first using a 7r/2 Bragg pulse 
to transfer half of the atoms from mode to mode 3 which 
is initially unpopulated. A weak Bragg pulse will not affect 
modes 1 or 2 due to the Doppler shift of light. 



B. Homodyne measurements 

The final step of the experiment is make a homodyne mea- 
surement on mode 1 using mode as the local oscillator, and 
similarly on mode 2 using mode 3 as the local oscillator. One 
can simultaneously apply two sets of Bragg pulses tuned to 
mix momenta fci with fco, and k2 with k^. 

To complete the homodyne measurements the population 
in each of the modes must be measured. The quadrature for 
mode 1 just before the Bragg pulse is proportional to the dif- 
ference in population found in modes 1 and after the pulse. 



{alit2)a,it2)V/^ 



{alit2)aoit2)V/' 



(12) 



where we have rescaled by the size of the local oscillator. Sim- 
ilarly for modes 2 and 3, 



X, 



a2{t2)al{t2)e^''^ + al{t2)a3{t2)e' 
{al{t2)a3it2)y/^ 

(4ft2)«3ft2))l/2 ' 



(13) 



As can be seen from the above equation, the ith quadrature 
angle depends on 9i, the relative phase between the relevant 
Bragg lasers. The phase is reversed in Eq. ( fT3] ) because k2 — k3 
is negative. The dependence on the relative phase has two im- 
portant consequences. Firstly, this parameter must be fixed 
from shot-to-shot in an experiment in order to measure the 
correct statistics of xf. Secondly, controlling 61 and 62 al- 
lows one to access quadratures of any angle in order to pro- 
duce a set of measurements that demonstrate entanglement or 
the EPR paradox. 



C. Entanglement criteria 

Measurements of the quadrature statistics can confirm the 
system is entangled by employing the appropriate separabil- 
ity or EPR criteria. In previous work ll4lll . we derived three 
entanglement criteria that take into account the quantum na- 
ture of the local oscillator. Earlier studies |42, 43, 44] were 
based on the simplified quadrature operators in Eq. (O, which 
in a quantum optics setting closely correspond to the mea- 
sured quadratures described by Eqs. ( fTSI ) and ( fT3] ). However, 
for the limited numbers of atoms employed in Bose-Einstein 
condensate experiments and the non-classical state of the local 
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oscillators generated by our proposed scheme, the difference 
between the simplified and measured quadrature operators is 
potentially important. The difference between earlier works 
and the criteria employed below is a direct consequence of the 
commutation relations of the measured quadrature operators, 





)=2^ 


(- 




(4«o> / 


^2 ' ^2 


)^2^ 


(- 


(4^2) \ 


{^1^3) J 



(14) 



where = XJ 

From this commutator, it follows that all separable states 
obey EH 



{bib,) 



2"2/ 



{blb^ 



(15) 



Violating this inequality proves the system is entangled. No 
state can violate the above inequality for all values of 4>i and 
02; generally these are experimentally adjusted to find the 
minimum value (i.e. maximal violation). We can express the 
left-hand side of Eq. ( fTSl l in the form 



LHS = 008^(0) (var[Xi - X^] + Var[Yi + K,] 

+ 4cos(<^)siii(0) {Na.r[X^,Y^] +Xa.v[Y^,X^] 
+ sin2(0) (Var[li + X^] + Var[Yi - ^2]) : (16) 



where the omitted superscript implies a quadrature angle of 
0, (01 +<^2)/2andVar[i,B] = \{AB) + \{BA) - 
{A) {13) is the covariance of A and B. Note that only the 
sum of the quadrature angles enters the expression; the value 
is independent of the difference 0i — 02- This simplifies the 
minimization problem to just one variable with a simple sinu- 
soidal form. 

For brevity, we define the separability parameter S, 



5 = 2 



Y&r[xt' ~Xp' 




l-{a\ix,)/{h\b,) 


+ 


1 - (4a2)/(&2^2> 



(17) 



The separability criteria is then simply 5 > 4. 

The Einstein-Podolsky-Rosen paradox is a stronger form of 
entanglement in the sense that it cannot be demonstrated by all 
mixed entangled states |49il. For our system, we previously 
showed ll4lll that violating 



Var[Xf ^ - X^^] + Var[y/^ + f/^] > 2 



1 - 



{a]a^) 



(18) 

for j = 1 or 2, demonstrates the EPR paradox. This is violated 
for at least one value of j if S* < 2. 



A similar, yet stronger EPR criterion is 

Yai[xf\X^'y 



Al,[X^^] =Var[x|^] 



ALf[^'2^^]=Var[f2*=] 



Vai[xf'] 
Var[Yf\Yi^f 



(19) 



Var[y/=] 
^ AUxt]AUY,^^] 

\-{a]a^)/{b]b^)'' ~ 

where in the last line we defined the EPR parameter E. Mini- 
mizing E with respect to 0i and 02 is less straightforward than 
for S, but the two will usually be minimized for similar phase 
angles'. 



IV. SIMULATIONS 
A. Outline of experiment 

We begin with a pure condensate of iVo = 10^ *^Rb 
atoms in a quasi- ID harmonic trap with trapping frequencies 
{ujzi tjj\^ = 27r X (1, 44) Hz. Due to the high trapping aspect 
ratio, we assume we can ignore the dynamics in the tightly 
trapped direction, and take a variational Gaussian ansatz for 
the transverse wave function. We have chosen system param- 
eters that closely match the previous work of Hilligs0e and 
M0lmer ll34ll to allow for a clear comparison of our results 
with their calculations. Such a quasi-lD condensate is de- 
scribed by the Lieb-Liniger Hamiltonian with an additional 
external potential. 



2m dz' 



+ V{z,t) + ^-y\z)^{z)^i;{z), 

(20) 

where m is the mass of an atom, g = 4:TTh'^as/mA± is an 
effective ID interaction constant, and V{z, t) = mcu'^z'^ /2 + 
Vl{z, t), where Vl is the optical lattice potential [Eq. ([T]i]. We 
use a scattering length Os = lOOao and obtain the effective 
cross-sectional area A± = 42 /im^ by minimizing the energy 
of a Gaussian ansatz, resulting in Il34ll 



A, 



2TTh 

miuj_ 



Vl + 2asniB, 



(21) 



where niD is the average linear density. The condensate 
ground state wave function is quasi-one-dimensional if niD S 
a^^. For these parameters the Thomas-Fermi length of the 
condensate is approximately 295 fim. 

The experimental sequence to generate and detect entangle- 
ment between two matter-wave pulses is illustrated in Fig. |2] 
To initiate the degenerate four-wave mixing, an optical lattice 



' Generally £'(</'i , <7i2) is minimal at the extreme points of S{(j>i, i}>2). where 
the (anti-)correlation is strongest. 
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the fragile entangled state. To avoid most of these collisions, 
the trapping potential may be ramped down adiabatically over 
a timescale short compared to the rest of the experiment (but 
fast compared to u;,-), removing most of the interaction energy 
from the system while preserving the quasi-one-dimensional 
state. This leaves the atoms to slowly expand radially in a 
Gaussian wavepacket. To model this procedure, in the simu- 
lations the nonlinear constant g is set to zero at time ti. 

The measurement scheme is implemented by applying a se- 
ries of Bragg pulses. The first pulse is also of wavelength 790 
nm, but is tuned to transfer half of the initial condensate at rest 
(in the lab frame) to mode 3 with momentum — 2fci and oc- 
curs between ti and t2- Immediately following this between 
t2 and ts are two simultaneous Bragg pulses of wavenumber 
(fco — fci)/2 = (fc2 — fc3)/2, each detuned to beam split modes 
1 with 0, and modes 2 with 3 respectively, as illustrated in 
Fig. 12(b). 



B. Simulation method 



FIG. 2: (Color online) (a) The intensity of the various optical lattices 
as a function of time. The solid (red) line indicates the intensity of 
the optical lattice generating the entanglement, which is ramped on 
and off with a continuous first derivative using parabolic curves. This 
process adiabatically transfers momentum to quasi-momentum states 
and vice-versa. The dashed (green) line indicates the Bragg pulse to 
split the phase reference in two. The dotted (blue) line indicates the 
simultaneous Bragg pulses used to beam-split each of modes 1 and 
2 with local oscillators in modes and 3. (b) A momentum space 
schematic of each stage of the process. 



of wavelength 790 nm and peak strength s = 1 is adiabatically 
ramped on and off over a time ti as depicted in Fig. EJa). In 
the simulations we ramp the lattice on and off using piece- 
wise parabolic curves with continuous first derivatives over 
a time period Tiamp = 2 ms. Gross-Pitaevskii simulations 
were performed to ensure that the shape and duration of the 
ramping curve efficiently transfers momentum states to quasi- 
momentum states and vice-versa. The optical lattice is de- 
tuned by 



(22) 



where ZVc = (43/64)fcL. In the frame of the lattice, the con- 
densate begins with mean momentum fco = ZVs, as depicted in 
Fig-Hi and degenerate four-wave-mixing generates new wave 
packets with mean momenta ki and fc2- 

At this point in the simulation collisions between atoms 
in the system become a hinderance. We therefore assume 
that the harmonic trap is removed after the optical lattice is 
switched off. During the subsequent expansion of the atom 
cloud the collision rate and hence the effective nonlinearity 
will decrease. Collisions have a negative impact on the mea- 
surement scheme and reduce the strength of correlation be- 
tween the quadratures. If the radial trapping is very strong, 
releasing the trap suddenly could potentially result in com- 
plex three-dimensional collision dynamics that could disturb 



We implement the truncated Wigner approximation ||38 
BTll to model the quantum dynamics of the proposed experi- 
ment. Despite being approximate, the method is numerically 
stable and has been shown to accurately treat quantum field 
dynamics on short to medium timescales. In particular, the 
validity condition of simulating a large number of particles 
compared to the number of modes is well satisfied in our nu- 
merical calculations lf52ll . 

The truncated Wigner method is implemented by stochas- 
tically sampling the initial Wigner functional W{ip{z)) and 
then evoling this according to the Gross-Pitaevskii equation 
(GPE) 



ih 



dipjz) 



2m dz^ 



+ V{z,t)+g\^{z)\^]^{z). (23) 



The initial state we use is a coherent state condensate found 
by solving the GPE in imaginary time with the addition of 
vacuum noise in the remaining empty modes ll52ll . The en- 
semble of trajectories then represents the evolution of the 
Wigner functional. Expectation values of symmetrically or- 
dered quantities are obtained by sampling moments of the 
field. For instance. 



i'*{z)^P{z) = {'4j\z)i>{z) + ^{z)^\z))/2. (24) 

We discretize the atomic field into 4096 points with a range of 
404 jiva, and the evolution is computed using a split-operator 
adaptive 9th order Runge-Kutta algorithm generated by the 
open-source software XMDS lUltl. We typically run 100-1000 
trajectories for each set of parameters. 



C. Results 

1. Maximal entanglement 

We begin with an analysis of simulations of the degener- 
ate four-wave mixing process only [up to time ti in Fig.|2a)] 
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FIG. 3: (Color online) Momentum density and relative number fluc- 
tuations versus the time held in the optical lattice, ti . (a) Population 
of each discrete momentum mode for the periodic model and (b) the 
trapped model. The initial condensate is barely visible in the peri- 
odic model as occupies just a single momentum mode. We see strong 
population growth around momenta fci and k2, as well as low energy 
scattering about ko. In (c) and (d), the solid line is the number in the 
entangled modes, {N-^ + iVj), and the dashed line is the variance of 
the number difference, Var[iVj — iVj]. Relative number squeezing 
is observed when this variance is less than the total number, and a 
strong correlation is exhibited at early times in (c), the data from the 
homogenous model, and (d), the data from the trapped model. At 
later times the non-classical correlation is lost. 



in order to identify the maximum violations of the entangle- 
ment inequalities without the complication of the measure- 
ment scheme. We employ both the trapped model described 
above and a simplified, periodic model, that does not include 
the axial harmonic trap. This allows us to illustrate important 
differences between the situations with well-defined quasi- 
momenta, and the more realistic trapped condensate. For the 
periodic model, we begin with 10^ atoms spread over 512 pe- 
riods of the optical lattice, or 202 /xm, which roughly corre- 
sponds to the high density region of the trapped condensate. 

In Fig |3ja,b) we plot the population of atoms in each 
momentum mode as a function of time held in the lat- 
tice, ti. This allows us to determine which momentum 
modes the degenerate four-wave-mixing process populates 
most strongly. The results show that the population grows 
fastest in the modes with momentum fco ± (105/256)A:l 
modulo 2kL for the periodic model, and ka ± (104/256)A:l 
modulo 2fci for the trapped model. These values are close 
to the predictions of the simple band-structure model, de- 
picted in Fig. [T] and differ because of the slight average 
density differences between the two models. For the re- 
mainder of this paper we use values of (fco, fci, fc2, ^3) = 



(0.672,0.262,-0.918, -1.381)fcL for the periodic model, 
and {ko,ki,k2,h) = (0.672, 0.266, -0.922, -1.381)fcL for 
the trapped model. 

The modes surrounding these peak growth modes also un- 
dergo significant growth. The time-energy uncertainty rela- 
tion allows momenta close to the resonance condition to ex- 
perience population growth for a timescale inversely related to 
their detuning. Below, we see that these additional populated 
modes do not prevent the detection of entanglement. After a 
longer period of time, secondary collisions populate a range 
of modes and the condensate is significantly depleted. Previ- 
ous experimental and theoretical work has shown this behav- 
ior is a result of the well-known d yna mical instability present 
in moving optical lattices ll36l[37[ Issll . Instabilities in optical 
lattices have traditionally been studied due to their thermaliz- 
ing effects; however the degenerate four-wave-mixing process 
we employ to create entanglement is a dynamical instability in 
mean-field terminology. 

In Ref. Isill . Hilligs0e and M0lmer implement a mean- 
field model and predict up to 95% transfer of the population 
to modes 1 and 2 (when an initial population seed is placed 
in mode 1). Our simulations account for spontaneous colli- 
sions into many modes and show that the transfer efficiency is 
actually limited to much less than this value without seeding. 
In this situation for the trapped case we reach a maximum of 
16% population transferred to the windows about fci and k2, 
and with a 10% seed we get a total of 45% conversion, show- 
ing that spontaneous scattering into multiple modes limits the 
transfer efficeincy. 

We have calculated the relative number squeezing between 
the two signal modes centered about momentum fci and fc2. 
In Fig.[3](c,d) we see the population in each mode grows ap- 
proximately exponentially as one would expect for this sys- 
tem ll45ll . It is important to note that a trapped condensate 
has a nonzero momentum spread. Detecting only a narrow 
momentum range of the trapped condensate would result in a 
poor overlap with the true spatial mode, degrading both the 
number correlation and the performance of the entanglement 
criteria. One must make measurements over a range of mo- 
menta of width 2dk. The number measured in signal mode j 
is then 

r+Sk 

TV. = / tjj^kj + k')'i}}{kj + fc') dk'. (25) 

J-Sk 

In Appendix lAl it is shown that the entanglement criteria de- 
scribed in Eqs. ( fT7lfT9] l hold true provided 5k is small enough, 
replacing (aja^) with {N^). 

In Fig.|4](c,d) we sum over five modes in the computational 
basis to determine the number of atoms in each pulse. We 
also see the number difference variance between the two sig- 
nal modes is significantly below the shot noise level {N^+N2) 
for smaller values of <i. The number squeezing degrades at 
later times due to secondary collisions transferring atoms out 
of the two signal pulses (in agreement with ll32ll ). 

We now analyze the entanglement of the modes 1 and 2 by 
implementing an idealized measurement scheme. In Fig. |4] 
we plot the results of direct measurement of the quadrature 
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FIG. 4: (Color online) Entanglement demonstration as a function 
of the hold time in the optical lattice, ti. The solid line represents 
the result for the homogenous calculation, while the dashed line is 
the result for the trapped calculation, (a) The separability criterion 
[Eq. J15H is maximally violated when the entangled populations 
A'^i + A^2 ~ 9. The presence of the trap reduces the level of vio- 
lation, but the system still demonstrates the EPR paradox as given by 
Eq. l |18l l. (b) The two systems display similar behavior to (a) with 
respect the EPR criterion, Eq. l |19t . Ensembles of 1000 stochastic 
trajectories were used to generate these results. 



TABLE I: The maximal amounts of separability and EPR violation 
predicted for periodic and trapped systems beginning with different 
initial seed populations. 



Number in Number 
Seed entangled difference Separability EPR 



Trap 


size 


modes 


variance 


violation 


violation 


No 





13 


0.5 


11.1 dB 


16.2 dB 


No 


10 


149 


16 


10.9 dB 


15.8 dB 


No 


10^ 


936 


124 


9.6 dB 


13.4 dB 


No 


10^ 


3920 


1053 


5.6 dB 


5.9 dB 


No 


10"* 


10465 


10243 


-0.1 dB 


-2.2 dB 


Yes 





9.3 


0.7 


4.5 dB 


3.5 dB 


Yes 


10 


38 


11 


4.5 dB 


3.5 dB 


Yes 


10^ 


296 


103 


4.4 dB 


3.3 dB 


Yes 


10^ 


2270 


986 


3.7 dB 


2.1 dB 


Yes 


10"* 


9963 


9836 


-0.1 dB 


-2.1 dB 



operators described by Eqs. (I12I13I I. The separability param- 
eter S and EPR parameter £ have been minimized with re- 
spect to the quadrature angles 0i and (f>2- To do this we di- 
rectly extract from the simulations the expectation values in 
the expression Eq. ( fT6l ). and from these determine the phase- 
angle (pi + (t)2 which minimizes S. We find that the optimal 
phase-angle varies significantly with time. For simplicity, we 
assume this phase-angle also minimizes £, which we find to 
be true in practice. Entanglement and the EPR paradox are 
indeed demonstrated for some time in both models. After a 
few milliseconds the quantum state of the field becomes com- 
plicated, and the entanglement detectable by these separabil- 
ity and EPR criteria disappears in agreement with previous 
work imiH. 

We include the magnitude of the maximal entanglement cri- 
teria violation (with respect to ti) in Table |T] along with the 



FIG. 5: (Color online) Minimal values of (a) the separability parame- 
ter S and (b) the EPR parameter £ with different initial seed popula- 
tions in mode 1 . Crosses (blue) represent results for periodic systems 
and circles (red) for trapped systems. Larger seeds appear to degrade 
the amount of entanglement violation. Each point corresponds to a 
different time ti that minimises the separability parameter for each 
system and seed size. 



number of entangled atoms. In both models the peak violation 
occurs when the total population in modes 1 and 2, N1+N2, is 
approximately 10. The turning points of S and £ correspond 
to when the number difference variance begins to grow. We 
note that the simplified, periodic system produces as much as 
10 dB stronger inequality violations compared to the trapped 
case. 



2. Seeding signal modes 

In optical experiments bright entangled sources can be cre- 
ated by 'seeding' the four-wave-mixing process with a co- 
herent input into one of the signal modes ifsill . We have in- 
vestigated populating mode 1 with a seed in order to boost 
the number of atoms in the entangled modes. This could be 
achieved experimentally with a short Bragg pulse before the 
optical lattice is applied. 

Figure|5]and Table |I]indicate the results of our analysis. In 
common with the optical case, we find that the entanglement 
violation is decreased for larger seeds. A seed of even 1% 
of the initial population, or 1000 atoms, results in noticeable 
degradation, and a 10% seed leads to no entanglement criteria 
violation at any time. We believe increased secondary scatter- 
ing is the cause of this degradation. Despite the fact that larger 
signals can be obtained by using a seed, the useable entangle- 
ment, such as measured by the entropy of entanglement lUTIl . 
does not increase with a larger seed^. 



^ This can be seen from Eq. (4); if one adds a coherent seed in mode 1 
by a displacement operation, d-^ —* Sj + a in the Heisenberg picture, 
then the resulting solutions are also simply displaced, and the entropy of 
entanglement is unchanged. 
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FIG. 6: (Color online) The measured dependence of the (a) sepa- 
rability and (b) EPR parameters on the phase of the Bragg lasers 
9i ^ 02 = <j}i + (t>2- Dots (black) indicate data from stochastic sim- 
ulations for each phase angle. Solid (red) lines represent fits to the 
data of (a) sinusoidal form with minimum value 5 — 1.52 and (b) 
two-frequency sinusoidal form with minima of £ = 0.43. These 
results are directly calculated from number correlations at time t^. 



3. Implementation of measurement scheme 

Finally, we have performed simulations implementing the 
sequence of Bragg pulses that are required to demonstrate en- 
tanglement in our scheme. To reiterate, first the optical lattice 
is applied to generate the entanglement until time ti. Immed- 
itately following this we set g and to zero (representing 
rapidly removing the trap on a timescale quick compared to 
1 /lOz but short compared to 1 /ur) and begin the Bragg pulse 
splitting the condensate into modes and 3. This pulse ends at 
time t2, when we apply two superimposed 7r/2 Bragg pulses 
which mix modes and 1, and modes 2 and 3. At time tj, 
the quadratures are measured by the appropriate number dif- 
ferences between the modes in fc-space. During the time the 
pulse sequence is applied the four separate condensates are 
overlapping in real space. To perform the homodyne mea- 
surements it would be necessary to allow these to expand and 
separate in real space, and then then accurately count the num- 
ber of atoms in each separate cloud. Note that a heterodyne 
scheme may be possible by investigating the real-space spec- 
tral components of the density at time ta, but may be difficult 
with limited detector resolution. 

We have carried out independent simulations of 360 dif- 
ferent Bragg pulse phase settings with 9i = —6*2. We find 
the optimal angle for violating the entanglement criteria in the 
same manner as would be necessary in an experiment. The 
results are plotted in Fig.|6] where each data point is the result 
of an ensemble of 100 stochastic trajectories. The statistical 
variations of these points give some indication of what may be 
expected from a similar number of experimental runs^. Previ- 
ous statistical experiments of BECs have involved tens lfl6ll to 
thousands ||7l!9[l of shots. 

We have performed a sinusoidal fit to the data to determine 
the maximal separability and EPR violation values. The sep- 



^ Although single and ensembles of truncated-Wigner trajectories are not 
formally equivalent to the expected results of experimental realizations, 
there is some degree of con'espondence Is^l . 



arability criteria follows a simple sinusoidal form [Eq. ( fTST ll 
with a minima of 1.52, or a 4.2 dB violation. The analytic 
form for the EPR parameter as a function of quadrature angle 
is more complicated; we fit the sum of two sinusoids (with pe- 
riods TT and 7r/2) resulting in the fitted minimum value of 0.43, 
or a violation of 3.8 dB. We see that these values fit closely to 
those of the idealized measurement scheme (c.f. Table|I]i, and 
that the measurement scheme as simulated is highly efficient. 
Of course technical issues such as detector efficiencies in a 
real experiment will potentially have a significant impact on 
the results. 



V. CONCLUSIONS 

We have combined the entanglement generating method 
of degenerate four-wave mixing for a trapped Bose-Einstein 
condensate loaded into a moving optical lattice proposed 
by Hilligs0e and M0lmer ||34|| with a quadrature measure- 
ment scheme and entanglement criteria previously derived in 
Ref. iiUl . Once the entanglement is generated, a series of 
Bragg pulses are applied in analogy to beam splitting used by 
homodyne schemes in quantum optics. Finally, the number of 
atoms in each mode is measured via time-of-flight expansion, 
and the quadrature values are determined from these. 

We have simulated a numerical model of a realistic quasi- 
one dimensional condensate that incorporates the effects of 
imperfect beam splitting by Bragg pulses, non-adiabaticity of 
the optical lattice, and the effects of the trapping potential. 
Our results indicate that neither multimode effects or the pres- 
ence of a trap will necessarily prevent the demonstration of 
entanglement with ultra-cold atoms. A theoretical analysis of 
different atom types, geometry and other experimental param- 
eters may lead to significant optimizations and improvements 
to the predicted entanglement criteria violations. 

These results are promising, but it is important to note 
potential experimental difficulties that we have neglected in 
our model. Most significantly, detector efficiency plays a 
crucial role in entanglement demonstration. Atom detection 
techniques are currently undergoing rapid development and 
we expect that modern detectors with efficiencies of 30-50% 
should be sufficient to demonstrate inseparability with this 
scheme llssll . It is critical that the phase of the Bragg pulse and 
optical lattice lasers be carefully controlled, in order to mea- 
sure the quadrature statistics of a specific quadrature angle. 
The effects of three-body loss and finite-temperature effects 
should also be taken into consideration. We have assumed 
the problem is purely one-dimensional; one may expect colli- 
sions into higher energy radial spatial modes during the four- 
wave mixing and the expansion process. These collisions can 
be reduced with tighter radial trapping frequencies, but truly 
ID systems suffer from phase fragmentation which may have 
other detrimental effects. 

Despite these potential complications we believe that we 
have outlined a feasible scheme for demonstrating entangle- 
ment between modes of ultra-cold bosonic atoms. We would 
like to point out that a similar detection scheme could be ap- 
plied to entanglement generated in a different manner, such as 
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by colliding two BECs in free space Il]J,|29l|30ll31l|32|]. How- 
ever, we remind the reader that the advantage of the quasi- 
ID geometry of the scheme we describe is it that limits the 
number of resonant modes and results in only two entangled 
matter-wave packets. It is possible to imagine other schemes 
that similarly populate few entangled modes; for instance a 
recent experiment observed well-separated and potentially en- 
tangled modes generated by four-wave mixing between dif- 
ferent hyperfine states of metastable helium jSoA . Not only is 
entanglement between massive particles interesting in its own 
right, but producing entangled atomic sources may lead to fur- 
ther novel experiments and quantum information applications 
in the future. 
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APPENDIX A: MULTIPLE MODES 



that the Bragg pulses act similarly for each momenta. We 
write the number measured about each relevant momenta 



r+Sk 

NAt)= '4>^ikj+k',t)dk', (Al) 

J-Sk 

and then measure each quadrature similarly to the single mode 
case. The numbers after the final Bragg pulse relate to the 
quadrature before; for example. 



-Sk 



^Hko + k',t2)ij^iki+k',t2)e-''^' 



-Sk 



+7/.t(fci + k\ t2)i^{ko + k', t2)e'^' dk' 

EE {N^{t2)Y'^Xt{t2). (A2) 

The above assumes 5k is within the range of values that 
the Bragg pulse affects, given by the time-energy uncertainty 
principle, h6k^/2m < 1/th- 

We can see that commutation relations are very similar to 
those in Eq. (fT4l) . 



It is straightforward to extend the entanglement criteria of 
Ref. iiHl to the situation where the signal modes have a finite 
momentum width. Realistic experiments will involve conden- 
sates with uncertain momentum, undergoing finite time-of- 
flight before measurement on detectors with non-zero reso- 
lution. 

We consider a range of momentum values small enough 



= 2i 1 - 



(^0(^2)), 



(A3) 



It follows that the entanglement criteria in Eqs. (I15I18I19I ) 
hold true, replacing (a]a^ ) with the total number in the mo- 
mentum range, (iV). 
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